Variable electrostatic transformer: controllable coupling of two charge qubits 
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We propose and investigate a novel metliod for tlie controlled coupling of two Josephson charge 
qubits by means of a variable electrostatic transformer. The value of the coupling capacitance is 
given by the discretized curvature of the lowest energy band of a Josephson junction, which can 
be positive, negative, or zero. We calculate the charging diagram of the two-qubit system that 
reflects the transition from positive to negative through vanishing coupling. We also discuss how to 
implement a phase gate making use of the controllable coupling. 



Following experimental demonstrations of individual 
qubits implemented with Josephson junctions operated 
in the charge [1,2] or flux [3-6] regimes, a lot of interest is 
now focused on building multi-qubit Josephson circuits. 
Recently, there has been an encouraging first experimen- 
tal step [7] demonstrating quantum-coherent dynamics 
of two charge qubits coupled directly through a capaci- 
tance. The coupling strength could not be varied in this 
experiment, since there is no simple way of changing the 
electrostatic capacitance between two metallic islands. A 
controlled coupling of charge qubits would be highly de- 
sirable. Many two-qubit transformations that could be 
implemented with a controllable coupling become more 
difficult or impossible at constant coupling, since in the 
presence of interactions one needs to correct for different 
dynamic phases of the two-qubit states. Two-qubit gates 
with untunable couplings have been suggested [8], how- 
ever, the gate design becomes significantly more complex. 

For flux qubits, the problem of a controlled coupling 
can be solved with variable flux transformers (see, e.g., 
[9]) which can be implemented naturally with the stan- 
dard tools of superconductor electronics [10]. For other 
types of "magnetic" qubits, e.g., spin qubits, the possi- 
bility to control the qubit coupling both in sign and ab- 
solute value has been predicted theoretically [11] (for re- 
cent experimental progress, see [12]). For charge qubits, 
however, a comparably natural solution has been lacking 
although the importance of the interaction control prob- 
lem has been recognized for some time. Suggested con- 
trolled coupling methods include direct spatial separation 
of qubit states by adiabatic transfer along junction ar- 
rays [13] or coupling through superconducting resonators 
[14-20]. The first one is too complicated to be practical 
at present. The second requires relatively large induc- 
tances of small geometric dimension that are fundamen- 
tally difficult to produce, in particular without creating 
additional dissipation at the relevant large frequencies. 

The purpose of this work is to analyze a new circuit to 
implement a controllable coupling of charge qubits. The 
circuit, shown in Fig. 1, generalizes the simple capacitive 
coupling used in [7] and makes use of the following basic 



feature of an individual small Josephson tunnel junction. 
If the junction dynamics is confined to the lowest energy 
band of its band structure, it behaves as a variable ca- 
pacitance with respect to the injected charge [21]. The 
magnitude of this capacitance depends on the bias point 
and can even be made negative. This is important for 
the application as a controllable coupling element, since 
it enables one to drive the coupling capacitance to zero. 
The physics of the capacitance modulation is the trans- 
formation of charge, when an injected charge causes a 
Cooper-pair transfer across the junction that changes the 
output charge. If, however, the two qubits are coupled di- 
rectly through a small junction (the coupling scheme con- 
sidered, e.g., in a somewhat different context in [22,23]) 
this mechanism would not produce a variable-capacitance 
coupling scheme, since variations of the junction capaci- 
tance are compensated by a redistribution of the trans- 
ferred charge on the capacitances in series with the junc- 
tion. Nevertheless, if the junction is included "perpendic- 
ular" to the coupling direction, as in Fig. 1, the required 
"variable electrostatic transformer" with a gate-voltage- 
controlled coupling capacitance is produced. 
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FIG. 1. Equivalent circuit of the variable electrostatic 
transformer for the controlled coupling of charge qubits. The 
Josephson coupling energy Ej of the small tunnel junction 
and the capacitances C, Cmi , Cm2 of the structure determine 
the dispersion relation of the junction. Its lowest energy band 
eo{q) provides the required variable coupling capacitance con- 
trolled by the gate voltage Vg. 

Indeed, if the voltage V is applied to the pair of in- 
put terminals (left or right, i = 1,2) of the circuit 
shown in Fig. 1, the Hamiltonian of the system is that of 
an individual Josephson junction (or "Cooper-pair box" 
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[24,25]) with the charge 2e{q — q') injected into it. where 
q = CmiV/2e and q' = CmiVg/2e. When the junction is 
in its lowest energy band eo, the output voltage Vout at 
the opposite terminals varies such that 



4eVCo = 2edVout/dq = dhoiq - q')/dq^ . 



(1) 



The input-output interaction strength defined by this re- 
lation varies with the charge 2eq' injected by the gate 
voltage and goes through zero at a bias point dependent 
on the shape of the energy band eo(g). 
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FIG. 2. Equivalent circuit of the two charge qubits coupled 
by a variable electrostatic transformer. 

If the transformer shown in Fig. 1 is inserted between 

two charge qubits (Fig. 2), it provides a gate- voltage- 
controlled qubit coupling. The qubits 1,2 with Josephson 
coupling energies Eji, i = 1,2, are biased by the gate 
voltages Vi through the capacitances Cgi and arc coupled 
to the central island of the transformer by capacitances 
Cmi- Inverting the capacitance matrix of this system, we 
can write the Hamiltonian as 

H = - Eji cos v-i + Eciirii - qi)^ 



1=1,2 



Ej cos If + Ec{n - q{ni,n2)f , 



(2) 



where Ec = 2e^l{C^ - Y^iClJC^i), Ea = 2e^/C^i, 
rii is the number of excess Cooper pairs on qubit i, 
and qi = ViCgi/2e. The qubits arc coupled through 
the charge q induced on the transformer junction: q = 
Qg-T,i{ni-qi)Cmi/CT:i, where qg = {Vg/2e) J2i Cmii'^- 
Cmi/CY,i)- The Josephson phase differences across the 
respective junctions are denoted by ipi and ip. Finally, 
C'si = Ci + Cgi + Cmi, and Cs = C + Cml + Cm2- 

We assume that the degrees of freedom of the trans- 
former junction are fast, and that the junction is confined 
to the lowest energy band of its band structure. Then we 
can replace the part of the Hamiltonian related to the 
transformer by the dispersion relation £0(9) of this low- 
est band, 

-E J cos ip + Ec{n - q{ni,n2)f ^ £0(5(^11, "-2)) • (3) 

This is the analogue of the Born-Oppenheimer approx- 
imation which in our case is valid if the characteristic 
energy gap between the bands of the transformer junc- 
tion is much larger than the qubit energies. For our 



qubit coupling, this condition requires that Ej ^ Eji 
for Ej < Ec, and {EcEjfl'^ > Eji for Ej > Ec (the 
qubits are assumed to be in the charging regime). 

To determine the qubit coupling provided by the trans- 
former we will use the known properties of the junction 
bandstructure [21,26]. Also, we use the fact that in the 
charge qubit regime {Eji <^ Eci, qi — 1/2), the dynam- 
ics of the charges rii in (2) is reduced to two states, so 
that the charges can be expressed through the Pauli ma- 
trices, rii = (czj + l)/2. Assuming that the structure in 
Fig. 2 is symmetric, i.e., Cmi/Csi = Cm2/CT.2 = c, the 
transformer term eo(g(ni,n2)) in the reduced Hamilto- 
nian can be expressed as follows: 

eo(g(ni,n2)) = V(7zi<yz2 + 5{u^i +0-^2) + A*- (4) 

The coupling coefficient v is given by 

^ = \ [^o{lo + c) + eo{qo - c) - 2eo(Q'o)] , (5) 

where qq = Qg + c^^{qi — 1/2). The term linear in 
shifts the qubit bias hy 6 — [£0(90 + 0) — €o('/o — c)]/4, and 
the constant term ^ = [eq (go — c) -I- €0 (go + c) + 2eo (go)] af- 
fects the energy of all two-qubit states and is not relevant 
as long as we are discussing one pair of qubits. 

Equation 5 is one of the central results of our work. 
We see that z/ has the structure of a discretized second 
derivative of £0 similarly to Eq. (1). Since £o(g) is a 
periodic function, u can be positive or negative. For 
instance, in the tight-binding limit Ec <C Ej, when 
eo(g) = — Acos(27rg), the coupling is: 



= — cos(27rgo)(l — cos(27rc)) , 



(6) 



and 6 = (A/2) cos(27rgo) cos(27rc). Obviously, the cou- 
pling can be controlled by the gate voltage through the 
average induced charge go, and can change sign. In gen- 
eral, the coupling constant can be calculated numeri- 
cally from Eq. (5). The results are plotted in Fig. 3 
which shows that in accordance with Eq. (6), the go- 
dependence of is harmonic for Ec < Ej, and changes 
sign at go ~ ±1/4. At larger Ec, J^(go) becomes non- 
harmonic, and for Ec S> Ej and small coupling coeffi- 
cients c, approaches c^Ec /2 for all go except for the vicin- 
ity of the point go = ±1/2, where v/c^Ec, becomes large 
in absolute value and negative: v/c^Ec — —Ec/2Ej for 



c < Ej/Ec, and vjc^Ec 



-l/2cfor EjjEc < c < 1. 



The periodic nature of £o(g) which has a minimum at 
g = and a maximum at g = ±1/2 implies that the max- 
imum of the absolute value of v is reached for c = 1 /2 and 
is equal to half the bandwidth, e.g., E'c'/S for Ec » Ej, 
and A for Ec<^Ej - see Eq. (6). 

The controlled electrostatic coupling of two qubits of- 
fers new possibilities to manipulate two-rjubit states. The 
simplest one is the change of the structure of the charg- 
ing diagram of the transformer-coupled qubits with the 



2 



gate voltage Vg . The charging diagram shows the regions 
of stability of the charge states (711,712) as a function 
of the induced charges g,, for small Josephson coupling 
Eji <IC Eci- The diagram is periodic with period 1 in 
both charges qi, and has a honeycomb structure in the 
case of fixed electrostatic coupling. Figure 4 shows ex- 
amples of the charging diagram. The diagram was ob- 
tained by minimizing the electrostatic qubit energy to- 
gether with the coupling energy (3). 
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FIG. 3. Coupling energy 1/ of the two qubits in units of the 

charging energy Ec of the transformer junction as a function 
of the average induced charge qo- The coupling strength is 
c = 0.1, and the ratio Ec/Ej is given for each curve. The 
normalization of v by (? makes it independent of c for small 
c. Inset: f(go) at qo = ±0.5, as a function of Ec/Ej for the 
same coupling strength c. 

At c = 0.3, the main qualitative features of Fig. 4 can 

be understood neglecting the difference between go and 
qg, from the go-dependence of the qubit interaction, see 
Fig. 3. For Qg = 0, the coupling is positive and the charg- 
ing diagram has the usual shape characteristic of a fixed 
electrostatic coupling. The cells of the honeycomb struc- 
ture are tilted in the qi = —q^ direction, since increasing 
the gate voltage of one qubit makes it easier to add a 
charge to the other qubit. In the absence of coupling, the 
degeneracy point of the four charge states (0,0), (0,1), 
(1, 0), (1, 1) is extended in the qi = q2 direction into the 
boundary between (0,1), (1,0), since positive coupling 
makes the charging energy of these two states lower than 
that of (0,0), (1,1). If the coupling is mostly negative, 
e.g., for qg = —0.5, the cells of the honeycomb struc- 
ture are tilted in the q\ = q2 direction. The degeneracy 
point of the four charge states is changed differently, since 
the states (0, 0), (1, 1) now have a lower charging energy 
than (0,1), (1,0). The degeneracy point is extended in 
the qi = —q2 direction and creates a direct boundary 
between (0,0) and (1,1). Around qg = 0.3, the qubit 
coupling is suppressed. Hence, four neighboring stability 
regions touch in one point, and their shape is approx- 
imately rectangular like for non-interacting qubits. In 
contrast to the case of fixed coupling, all the cells in Fig. 4 



arc curved, since changing qi, 52 also changes the effective 
charge qo of the transformer island, so that the effective 
qubit coupling constant v and the renormalization 5 of 
the qubit bias change with qi. This effect is largest for 
values of qg around which v{ciq) varies strongly, like e.g., 
qg = —0.4, see Figs. 3 and 4. 




FIG. 4. Charging diagram of two charge qubits coupled 
through the variable electrostatic transformer for c = 0.3 and 
Ec /Ej = 3. The curves show the boundaries of the stability 
regions of the charge states (711,712) in the plane of the two 
induced charges qi = ViCgi/2e, i = 1,2. The gate charge qg is 
controlled by the transformer gate voltage Vg. The changing 
shape of the stability region reflects the Vg-controUed transi- 
tion from positive to negative effective coupling capacitance 
of the two qubits. Dashed lines: uncoupled case. 

The possibility to switch the coupling between two 
qubits on and off is important for the realization of prac- 
tically all types of two-qubit gates. An advantage of 
our scheme is that the effective qubit coupling can be 
suppressed even in presence of a small direct positive 
geometric capacitance between the qubit islands. This 
"parasitic" capacitance only shifts the values of the gate 
voltages at which the total coupling vanishes. As an ex- 
ample, we show how our proposal can be used to directly 
realize a phase gate, a transformation that changes the 
sign of the |11) state and does not change the amplitudes 
of the other three two-qubit states. If the qubit tunnel- 
ing amplitudes are tuned to zero, the Hamiltonian of the 
two-qubit system (Fig. 2) with coupling (4) is 

F = jycr^l(T22 + 7/(o-zl + ^^2) , (7) 

where 77 = 6—{qi — l/2)Eci- (For simplicity, we again dis- 
cuss only the symmetric structure, which for the purpose 
of realizing the phase gate also implies identical qubit 
charges qi. All results of our work can be directly ex- 
tended to the asymmetric case.) By adjusting the three 
gate voltages Vi, V2, and Vg such that u{qo = q*) = 
and qi = 1/2 + 5{q*)/ Eci one can make the Hamiltonian 
(7) vanish. For Ej <C Ec such a stationary qubit state is 
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realized when qi = 1/2 and qq = Qg = 1/4. By applying 
gate-voltage pulses one can temporarily switch on both 
the interaction and the qubit bias in such a way that the 
qubit states accumulate dynamic phases [27]. To obtain 
a phase gate, i.e., to have the state |11) accumulate the 
phase TT up to an irrelevant common phase of all four 
qubit states, the phase (j)^ = J dtv{t) due to the interac- 
tion energy and the phase <Pti ~ I dtri{t) due to the qubit 
bias should be chosen as 4>v = 4'r)= 7r/4. 

These conditions fix only the time integral of the gate- 
voltage pulses. Their shape can be chosen to reduce 
the excitation amplitude of the upper energy band of 
the transformer junction that would invalidate our Born- 
Oppenheimer approximation (3) and effectively entangle 
the qubit states with the states of the transformer creat- 
ing an additional decoherence mechanism. Transitions to 
the upper energy bands are suppressed if the pulses are 
sufficiently slow on the scale of the transformer energies. 
In this case, the amplitude ak (t) of an excitation process 
to the kth band can be calculated by standard adiabatic 
perturbation theory as 

a.= r..i^^^^e-/:-(^')-', (8) 
J ek{T) 

where is the energy of the kth band relative to cq. 
Since the total phases accumulated by the qubit states 
during the gate-voltagc pulses arc of the order of tt, the 
adiabatic condition also means that the pulse amplitude 
should be small, i.e., qo = q*+x{t), a; ^ 1. Therefore, the 
energies and matrix elements in Eq. (8) can be evaluated 
at qo = q* , and Eq. (8) reduces to 

|afc|=2£c|(fc|(n-g)|0)||x(efe)|. (9) 

Here, x(efe) ~ J dTx{T)e^'^'''^ is the Fourier component 
of x{t) which decreases exponentially in CfeTo for smooth 
pulses, where ro is the characteristic pulse time. For 
Ej Ec, the energy-band gap of the transformer is of 
the order of Ej, and this sets a lower limit for the pulse 
time through Eq. (9). For Ej > Ec, the gap between 
the energy bands grows and the adiabaticity condition is 
easier to satisfy (see the discussion after Eq. (3). 

In conclusion, we have proposed a method to couple 
two Josephson charge qubits by an electrostatic trans- 
former that produces an effective coupling that can be 
varied in magnitude and sign. We have given an explicit 
expression of this coupling as a discretized second deriva- 
tive of the energy band of a small Josephson junction. 
Our proposal works even for asymmetric structures and 
in the presence of parasitic capacitances. It allows the 
implementation of a variety of two-qubit gates, and we 
have explicitly demonstrated how to build a phase gate. 
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